Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN ETjiX style file v2.2) 



Systematics in the X-ray Cluster Mass Estimators 



E. Rasia 1 , S. Ettori 2 , L. Moscardini 3 , P. Mazzotta 4 5 , S. Borgani 6 ' 7 8 , K. Dolag 9 , 
G. Tormen 1 , L.M. Cheng 10 , A. Diaferio 11 



O 
O 

(N 



1 Dipartimento di Astronomia, University di Padova, vicolo dell' Osservatorio 2, 1-35122 Padova, Italy (rasia,tormen@pd. astro. it) 

2 INAF, Osservatorio Astronomico di Bologna, via Ranzani 1, 1-40127 Bologna, Italy (stefano.ettori@bo.astro.it) 

3 Dipartimento di Astronomia, Universita di Bologna, via Ranzani 1, 1-40127 Bologna, Italy (lauro.moscardini@unibo.it) 

4 Dipartimento di Fisica, Universita di Roma Tor Vergata, via delta Ricerca Scientifica 1, 1-00133 Roma, Italy (pasquale.mazzotta@roma2.infn.it) 
, 5 Harvard-Smithsonian Centre for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 

■ 6 Dipartimento di Astronomia, Universita di Trieste, via Tiepolo 11, 1-34131 Trieste, Italy (borgani@ts.astro.it) 
JjT^ 7 INFN - National Institute for Nuclear Physics, Trieste, Italy 

8 INAF — Osservatorio Astronomico di Trieste, via Tiepolo 11, 1-34131, Trieste, Italy 

9 Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild Strasse I, D-85748 Garching bei Miinchen, Germany (kdolag@mpa-garching.mpg.de) 

10 Institute of Theoretical Physics, Chinese Accademy of Sciences, Beijing 100080, China (clm@itp.ac.cn) 

11 Dipartimento di Fisica Generate "Amedeo Avogadro", Universita di Torino, via Giuria 1, 1-10125, Torino, Italy (diaferio@ph.unito.it) 



o 

(N 



O 

CO 



> 

1st draft: jul 05 

<N 

O ' ABSTRACT 

We examine the systematics affecting the X-ray mass estimators applied to a set of five 
^ (— i galaxy clusters resolved at high resolution in hydrodynamic simulations, including cooling, 

O star formation and feedback processes. These simulated objects are processed through the 

X-ray Map Simulator, X-MAS , to provide Chandra -like long exposures that are analyzed 
to reconstruct the gas temperature, density, and mass profiles used as input. These clusters 
have different dynamic state: we consider an hot cluster with temperature T = 11.4 keV, a 
perturbed cluster with T = 3.9 keV, a merging object with T = 3.6 keV, and two relaxed 
systems with T = 3.3 keV and T = 2.7 keV, respectively. These systems are located at 
z = 0.175 so that their emission fits within the Chandra ACIS-S3 chip between 0.6 and 1.2 

-^500- 

We find that the mass profile obtained via a direct application of the hydrostatic equilib- 
rium equation is dependent upon the measured temperature profile. An irregular radial dis- 
tribution of the temperature values, with associated large errors, induces a significant scatter 
on the reconstructed mass measurements. At -R2500, the actual mass is recovered within 1 a, 
although we notice this estimator shows high statistical errors due to high level of Chandra 
background. Instead, the poorness of the (3— model in describing the gas density profile makes 
the evaluated masses to be underestimated by ~ 40 per cent with respect to the true mass, both 
with an isothermal and a polytropic temperature profile. We also test ways to recov er the mass 
by adopting an an alytic mass model, such as those proposed by Nav arro et alJ i ll 9971) and 
Rasi a et alT l l2004l) . and fitting the temperature profile expected from the hydrostatic equilib- 
rium equation to the observed one. We conclude that the methods of the hydrostatic equilib- 
rium equation and those of the analytic fits provide a more robust mass estimation than the 
ones based on the f3— model. In the present work the main limitation for a precise mass re- 
construction is to ascribe to the relatively high level of the background chosen to reproduce 
the Chandra one. After artificially reducing the total background by a factor of 100, we find 
that the estimated mass significantly underestimates the true mass profiles. This is manly due 
(i) to the neglected contribution of the gas bulk motions to the total energy budget and (ii) to 
the bias towards lower values of the X-ray temperature measurements because of the complex 
thermal structure of the emitting plasma. 

Key words: cosmology: miscellaneous - methods: numerical - galaxies: clusters: general - 
X-ray: galaxies - hydrodynamics. 
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1 INTRODUCTION 

The gravitating mass in galaxy clusters is the most fundamen- 
tal quantity to use them as cosmological probes. Measurements 
of the total cluster masses have been used for the last 70 years 
(since lZwickvl 1933b to infer the presence of dark matter over Mpc 
scales and to constrain the cosmological parameters that describe 
the distribution and growth of the cosmic density fluctuations of 
> 10 14 M(T) viri alized structures (see the reviews by iRosati et alJ 
l2002tlVoItl 2005). These structures appear as well defined and eas- 
ily observable at the X-ray wavelengths, where the emission is 
mainly dominated by bremsstrahlung processes and is proportional 
to the square of the plasma density. This observational evidence al- 
lows to identify massive relaxed clusters and to study their physical 
properties through the estimates of the X-ray emitting intra-cluster 
medium (ICM) density and temperature. Since the sound crossing 
time in the ICM is sufficiently shorter than the age of the structure, 
the assumption that the ICM is in hydrostatic equilibrium within the 
underlying dark matter potential generally holds (Sarazin 2002). 
ICM density and temperature are used to determine the radial mass 
profile via the hydrostatic equation (see, e.g.. lCowie et aljfl987h . 
with a further assumption of the spherical symmetry for the sake of 
simplicity. 

In the present work, we investigate the systematics that affect 
the X-ray mass measurements by using hydrodynamic simulations 
proce ssed by our X-ray Map Simulator (X-MAS , iGardini et alJ 
2004) to produce realistic event files and images that are then ana- 
lyzed in a way identical to what is done with observed clusters sys- 
tems. Several authors studied the robustness of the hydrostatic equi- 
librium and the uncertainties related to this mass estimator. General 
conclusions were that there is a good agreement between X-ray 
and true masses, in particular in the outer regions, wi th a stan- 
dard deviation of approximately 15 per cent (see, e.g., ISchindleJ 
199dlEvrard et al J 19961: cfr . also lBartelmann & Steinmetzll996l) . 
Balland & Blanchar f F997l) )ointed out, however, that the uncer- 
tainties on the mass measurement propagated from the weak con- 
straints on the gas temperature profile via the hydrostatic equi- 
librium equation make any accuracy on the determination of the 
mass profile very poor. More recently , iRasia et alJ 12004T) (see also 
iKav et aljEool iBorgani et ail [2004) found in a set of hydro-N- 
body simulations that the ICM is not in a perfect hydrostatic state 
and, thus, the masses can be underestimated by up to 20 per cent 
in non-radiative ICM models due to residual gas bulk motions. 
IKav et all 120041) conclude that more thermalization (and less sig- 
nificant departures from the true mass estimates) is obtained when 
cooling and feedback are included and that an isothermal /3-model 
still provides accurate mass estimates within 20-30 per cent un- 
certainties. With our approach in which (1) the known input on 
the matter distribution is properly convolved with the response of 
real instruments and (2) an observational-like analysis is performed 
on the output to recover the original input, we extend the previous 
work limited to the study of the hydrodynamic simulations in or- 
der to highlight the systematic effects present both in the assump- 
tion of hydrostatic equilibrium and in the observational techniques 
adopted in recovering the mass estimates. 

The paper is organized as follows. We introduce our simulated 
dataset in Section 2. The outputs produced by combing the hydro- 
dynamic simulations and the X-ray Map Simulator for Chandra 
data are presented in Section 3 and analyzed to recover the gravi- 
tating mass in Section 4. The results are presented in Section 5 and 
summarized and discussed in Section 6. All the errors quoted are 



at la level (68.3 per cent level of confidence for one interesting 
parameter). 



2 THE SAMPLE OF SIMULATED CLUSTERS 

The simulated clusters have been taken from two different simula- 
tions in order to obtain a sample with a significant range of masses 
and dynamical states. 

The simu lations has been carried out with GADGET-2 
lSpringell2005l). a new version of th e parallel Tree-SPH simulation 
code GADGET (Snringel et al. 2001). It uses an en tropy-conserving 
formulation of SPH ( Soringel & Hernauist 2002), and includes ra- 
diative cooling, heating by a UV background, and a treatment of 
star formation and feedback processes. The latter is based on a sub- 
resolution mode l for the multi-phase structure of the interstellar 
medium I S pringel & Hernquistl2003t) . 

The first cluster, Cnot, has been extracted from a dark-matter 
only simulation with a box-size of 479 hr 1 Mpc lYoshida et alJ 
2001). The assumed cosmological model is a standard flat ACDM 
universe, with f2 m = 1 - f2 A = 0.3, a$ = 0.9, Q, h h 2 = 0.019 
and H = lOO/i km s _1 Mpc -1 with h = 0.7. This cluster has 
been resimulated at higher mass and force resolution. The new ini- 
tial conditions for this system have been generated by applying 
the Zoomed Initial Condition (ZIC) technique iTormen et all 19971) . 
This method allows one to increase the mass resolution in a suitably 
chosen high-resolution Lagrangian region surrounding the struc- 
ture to be re-simulated, and at the same time to correctly describe 
the large-scale tidal field of the cosmological environment by using 
low-resolution particles. The mass resolution of the gas particle for 
this cluster is 1.7 x 10 s M@, the gravitational softening length 
was kept fixed at e = 30.0/i _1 kpc comoving (Plummer-equivalent) 
and was switched to a physical softening length of e — 5.0ft _1 kpc 
at z = 5. The resimulation foll ows star formation, feedback and 
heatin g by thermal conduction Jjubelgas et aljEo04l iDolag et"all 
120041) . The SN efficiency in powering galactic winds is set to 50 
per cent which turns into a wind speed of 340 km s _1 . 

The other four simulated systems, Cp er t, C'M erg , Cn e n, and 
CR e i2, used in the following analysis, have been extracted from the 
large- scale hydrodynamic simulation described in Borga nTet alJ 
(2004). We refer to that paper for a detailed description of that 
simulation, while we provide here only a short summary. The sim- 
ulation follows the evolution of 480 3 dark matter (DM) particles 
and an initially equal number of gas particles within a box of 
192 /t -1 Mpc on a side of the previous cosmological model, but 
with a lower normalization of the power spectrum, as = 0.8. 
The mass of the DM and gas particles are, respectively, tti-dm = 
4.6 x 10 9 /i _1 M Q andm gas = 6.9 x 10 8 /i _1 M Q . The Plummer- 
equivalent gravitational softening was kept fixed at epi = 7.5 h~ l 
kpc comoving and switched in physical units at z = 2. 

Cpert is directly extracted from the final output of this sim- 
ulation by using the standard identification criterion based on the 
spherical over-density. From the same simulation we selected the 
remaining three objects which, in an analogous way of Cnot, have 
been re-simulating at high-resolution, following the ZIC tecnique. 

These runs, wh ose results are extensively discussed elsewhere 
( Borgani et al. 2005), have been performed by assuming a mass for 
the DM and gas particles 10 times smaller than in the original sim- 
ulation; also the Plummer-equivalent gravitational softening has 
been reduced and corresponds to epi = 3.5 h^ 1 kpc at z = 0. For 
these new runs the feedback scheme is calibrated to lead to higher 
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Table 1. Properties of the simulated galaxy clusters forming our dataset 
(their identification names are given in the corresponding columns). i? v j r 
is the virial radius; R500 and -R2500 are the radii within which the over- 
density with respect to the critical value is 500 and 2500, respectively; M vlT 
is the virial mass. JVdm and N gaE represent the number of DM and gas 
particles inside B v ir> respectively. T s i is the spectroscopic-like temperature 
inside i?, v j r , calculated by summing over the particles having T > 0.5 keV 
(T > 1 keV for C Pert ); r SiNPW and c NFW and r SiRTM and c RTM are 
the scale radius (rescaled at the observation redshift) and the concentration 
of the NFW and RTM mass model, respectively. 





CHot 


Cpert 


CMerg 


Cpe.ll 


CrcU 


Rvir [kpc] 


2713 


1967 


1662 


1567 


1368 


R500 [kpc] 


1255 


1225 


808 


732 


646 


-R2500 [kpc] 


515 


585 


361 


305 


265 


M vir [1O 14 M ] 


22.6 


7.0 


4.1 


3.6 


2.3 


AT DM [10 5 ] 


17.4 


0.9 


5.5 


4.6 


3.0 


Afgas [10 5 ] 


14.1 


0.7 


4.2 


3.9 


2.6 


T sl [keV] 


11.4 


3.9 


3.6 


3.3 


2.7 


l"j,NFW 


438 


250 


182 


168 


142 


CNFW 


4.5 


5.0 


5.3 


5.4 


5.3 


fs,HTM 


170 


93 


63 


57 


46 


c rtm 


11.7 


13.5 


15.6 


16.4 


16.4 



galactic wind velocities with respect to the original feorgani et alJ 
bOOd ' simulation (« 480 km/s instead of » 360 km/s). 

The main properties of the five clusters forming our dataset 
are summarized in Table Q In particular we list: the virial radius, 
i? v i r , defined as the radius at which the over-density assumes the 
value dictated by the spherical top-hat model (see, e.g., lEke et alJ 
119961) : the radii -R2500 and Rsoo, within which the over-density 
with respect to the critical value is 2500 and 500, respectively; 
the virial mass M vlr , i.e. the mass included inside i? v j r ; the num- 
ber of DM and gas particles inside the virial radius (TYdm and 
iVgas, respectively); the spectroscopic-like temperature, defined 
as T s i = J WTdV/ J WdV, where the weighting function W 
equals to n 2 /T ' 75 , being n and T the density and temperature of 
each gas particle. This temperature definition has been introduced 
in the ana lysis of the hydrodynamic simulations of galaxy clus- 
ters by Mazzotta et al. 1 2004) to provide a better approximation to 
the values extracted from fits of observed X-ray spectra (see also 
Vikhli nirj2005l for an extension to cooler systems considering the 
effects of metals). 

To summarize, the following five objects have been selected 
as examples of clusters with different thermal and dynamic states: 

• C Ho t (T s i = 11.4 keV, M vlr = 2.26 x 10 15 M©): a forming 
cluster, still accreting mass from the outskirts; 

• Cpert (T s i = 3.9 keV, M vir = 7.0 x 10 14 M Q ): a perturbed 
cluster which shows in the temperature map a cold substructure 
infalling toward the centre; 

• C M erg (T s i = 3.6 keV, M v ir = 4.1 x 10 14 M Q ): an object 
that experienced a recent major merger; 

• Cfleii (T s i = 3.3 keV, M vil = 3.6 x 10 14 M ): a relaxed 
structure; 

• C Re i2 (T s i = 2.7 keV, M vir = 2.3 x 10 14 M Q ): a colder 
relaxed cluster. 



3 MOCK X-RAY OBSERVATIONS OF THE 

HIGH-RESOLUTION SIMULATED CLUSTERS 

The sample of our hydro-N-body simulated galaxy clusters has 
been processed with X-MAS to obtain mock Chandra ACIS-S3 ob- 
servations. The event files produced with our program are then an- 
alyzed in the same way and with the same tools of real observa- 
tions. A detailed description of X-MAS is reported in Gardi nTet alJ 
(2004). The design of the two main units constituting the software 
is here summarized. 

• First Unit: generation of differential flux maps. The input 
of this unit is the output of a hydro-N-body simulation. For each 
gas particle, we compute the emissivity and distribute it over the 
corresponding volume. After selecting a line of sight for the sim- 
ulated observation, we compute the projected spectrum for each 
pixel and then the differential flux for each angular coordinate in 
bins of energy. The final step is to add Galactic absorption. 

• Second Unit: simulation of Chandra observations. We 
estimate the expected number counts and iteratively subdivide 
the tile region until the counts become smaller than a given 
threshold (10 counts); we then use the command FAKEIT in 
the utility XSPEC (see, e.g., Xspec User's Guide version 12. 2 1 ; 
iDorman et al Jl2003h to convolve the spectral model of each sub- 
region with the response of the CCD and to add the sky back- 
ground (acis_c_S3_bg_evt_191000.fits); at the end we generate the 
final photon event file. 

A Galactic equivalent column density of Np — 5 x 10 20 
cm -2 is assumed and the cluster metallicit y has been fixed at 
Z = .3 Zq, using the solar value tabulated in lAnders & Grevessd 
(1989). In order to be able to observe a physical size of ~ 1.5 Mpc 
within the field of view of Chandra ACIS-S3 (8.3 arcmin) we im- 
pose a redshift z = 0.175 for the simulated clusters. 

Since we want to study the systematic discrepancies between 
observed and real quantities for ideal observations, we have ap- 
plied very long exposure times to all the simulated observations 
aiming to minimize the statistical uncertainties related to the num- 
ber counts: 0.5 Msec for Cpot, 1 Msec for Cp er t, Cuierg, Cne.ii 
and Cr b i2- Note that, despite the high exposure time, the effective 
net counts number is of order of 20 per cent or less in the bins more 
external than -R2500 due to the fact of the background. 

In the left panels of Fig. Q we report the photon images of 
our simulated sample. All the images are binned to 2" pixels and 
divided by their instrument maps. They are computed in the soft 
[0.5-2 keV] energy band, to better resolve the presence of small 
high-emissivity structures or merging objects that otherwise could 
be hidden by the background. The central panels of Fig.Qpresent 
the spectroscopic-like temperature maps, with the isoflux contours 
over-plotted. These maps are directly produced from the gas parti- 
cles in the simulation and show the thermal structure of the systems 
as would been obtained by ideal X-ray observations of the simu- 
lated clusters. On the right panels of the same figures, we plot the 
histograms of the emission measure, i.e. the sum of the square of 
the gas density of all the particles inside the box image for Chou 
Cp e rt and C M erg or inside -R2500 for C Re ii and C Re i2, for each 
bin of temperature and per unit volume. The vertical lines repre- 
sent the cluster temperature, Tisoo- We notice that the distribution 
is quite regular and Gaussian only for the relaxed cluster C Re i2, 
while many features are evident for the other systems, revealing 
their perturbed dynamic state. Nevertheless, we notice that T2500 is 

1 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/XspecManual.pdf 
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Figure 1. Images obtained from the simulated galaxy clusters: different rows refer to C/fot. Cperi. CWera! Crm an d Crc12> from top to bottom. Left 
column: photon images in the X-ray soft [0.5-2 keV] energy band. The images are 8.3 arcmin-width, exposure-corrected, and binned to 2"-pixel. The green 
circles show the regions masked out in the analysis. Central column: X-ray logarithmic isoflux contours over-plotted to the spectroscopic-like temperature 
map, both extracted directly from the hydrodynamic simulations; the temperature scale is shown on the left. Right column: emission measure from the region 
inside R500 (inside the box photon image for Cnot, Cp er t); the solid vertical black line refers to T2500 (see Tablel2l. 
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always lower than the value corresponding to the peak in the emis- 
sion measure distribution. This is the consequence of fitting with 
a single temperature the spectra of cluster regions that, as high- 
lighted by the temperature map, s how significant gradie nts or are 
highly thermally complex (see, e.g. lMazzotta et alj2004) . 



3.1 Description of the clusters in the sample 

Cfioi and Cp e rt have a quite complex cluster structure. This is ev- 
ident in the emission measure distributions, which are very broad. 
The fact that they extend to low temperatures indicates the presence 
of many cold blobs of gas. Some of them can be seen in the soft X- 
ray images as bright X-ray blobs distributed around the clusters. 
The cool nature of these X-ray blobs is confirmed by the tempera- 
ture maps. To estimate the cluster mass of these complex clusters, 
we follow the observational approach of excluding from the spa- 
tial and spectral analysis all sub-clumps clearly detectable from the 
soft X-ray image. The excluded regions are shown as green circles 
in Fig.Q We know that this procedure partially helps to reduce ther- 
mal contamination from colder substructures, however many ther- 
mal structures may still remain as they do not produce detectable 
perturbations in the soft X-ray image. 

The largest and hottest cluster, Chou is actually a still forming 
object that has not reached the equilibrium state and is undergoing 
merging events by small substructures. 

Cpert presents, in the southern hemisphere, a sub-clump in- 
falling toward the cluster centre and merging with a cold structure 
located at 1.5 arcmin from the centre. 

Cnierg appears very relaxed, with no evidence of any peculiar 
structure both in the photon image and in the temperature map. The 
elliptical shape of the isoflux contours is the only imprint of a recent 
major merging happened at z ~ 0.2. A further evidence of this 
recent impact with another big clump comes from the histogram 
of the emission measure that peaks at two different temperatures 
(T ~ 4.1 keV and T ~ 3.4 keV). 

The simulated clusters Cpeii and Cp e i2 have a regular photon 
image with round isophotes. On the other hand, the temperature 
map of CR e i2 shows the presence of a cold arc in the northern part 
at w 1 arcmin from the centre. 



4 AN OBSERVATIONAL-LIKE X-RAY ANALYSIS 

The event files produced by X-MAS have been then processed using 
the CIAO 3.0.1 software (see http://cxc.harvard.edu/ciao/ 1 and the 
calibration files in CALDB 2.23. Surface brightness profiles and 
spectra have been created and analyzed as described below. In both 
analyses we exclude the sub-clumps detected in the photon image 
of Cnot and Cp er t- Moreover, in all clusters we mask the inner 50 
kpc to avoid the influence of the central cooling part. The excluded 
regions are marked by green circles in the photon images. 



4.1 Spatial analysis: brightness profiles 

The surface brightness profiles have been extracted from [0.5-5 
keV] images that have been corrected thanks to the corresponding 
instrumental maps. 

To build the profile we consider annuli centred on the X-ray 
peak (in any case, we verified that the off-set from the actual min- 
imum of the cluster potential well is never larger than 2 arcsec) by 
requiring at least 5000 total counts per bin. 



Table 2. Best-fit parameters used for the X-ray mass estimators and ob- 
tained from the spatial and spectral analysis of the five objects in our sam- 
ple. r c and /3 are the core radius (in kpc) and the power coefficient of the 
/3— model; cj,kg is the constant representing the background (Eq.0; 7 is the 
polytropic index; X 2 e< i ^ and X 2 e< i are the reduced x 2 associated to the 
/3-model and polytropic relation fitting, respectively; d.o.f. are the degrees 
of freedom of the fits; T2500 is the spectroscopic temperature calculated 
inside R2500; ""2500 is the associated error at 68 per cent confidence level; 
X^ ed is the reduced x 2 , evaluating the goodness of the spectral fit; r a NFW 
and cnfw aR d r K rtm an d cjtm are the scale radius (in kpc) and the 
concentration of the NFW and RTM mass model, respectively. 





Cnot 


Cp e rt 










85±2 


101±6 


61±2 


21=1=1 


14=1=1.5 


P 


0.56±0.01 


0.44±0.01 


0.62±0.01 


0.52±0.01 


0.56±0.01 


c bkg 


1.78 


4.30 


3.91 


3.88 


4.06 




4.2 (13) 


3.9 (70) 


7.6 (44) 


5.3 (39) 


11.7 (29) 


7 


1.09±0.03 


U1±Q.15 


1.10±0.02 


1.11 ±0.04 


1.06±0.03 


X 2 j (d.o.f.) 


1.4 (6) 


0.6 (4) 


0.5 (10) 


0.9 (5) 


2.6(3) 


^2500 


11.18 


3.87 


3.76 


3.39 


2.64 


CT2500 


±0.48 


±0.36 


±0.08 


±0.12 


±0.08 


X 2 .^ (d.o.f.) 


0.51 (45) 


0.36 (39) 


0.44 (103) 


0.31 (30) 


0.37 (80) 


r s ,NFW 


668±47 


350±100 


216±103 


82±45 


252±38 


cnfw 


6.0±0.3 


6.1±2.5 


8.6±2.3 


16.2±7.9 


7.2±1.0 


r s,RTM 


693±112 


351=1=121 


148±69 


41 ±20 


226±70 


CRTM 


7.5±1.1 


7.4±5.3 


14.3±5.4 


36.2±8 


9.9±3.3 



The azimuthal profile has been fitted with a /3-model 
ICavaliere & Fusco-Femianoll976t) : 



S(r) = So 



1 + 



(0 



0.5-3/3 



' c bkg, 



(l) 



where r c is the core radius, f3 is the power coefficient, and the con- 
stant Cbkg is the value added to take into account the background 
present in the image. The best-fit parameters obtained by analyzing 
our clusters are listed in Table|2| For all our systems, the deviations 
between the surface brightness profile and the /3-model fit are less 
than 10 per cent in the interval between 0.1i? v ir and li? v ir (see, 
e.g., Fig.|2|for cluster Cp e n). Nevertheless, the model shows a sig- 
nificant discrepancy in unit of a, as shown by the large values of 
the reduced \ 2 reported in Table [2] Overall, we conclude that the 
j3— model cannot reproduce the profile of the X-ray emission of 
the simulated clusters and we find that also a double /3-model is 
statistically excluded by our data. 



4.2 Spectral analysis 

To measure the overall temperature we extract spectra from circular 
regions centred on the cluster centre and with a radius r = -R2500 ■ 
To calculate the temperature profiles we extract spectra from annu- 
lar regions for which the net counts (i.e. belonging to the cluster) 
are greater than 5000 and are at least 25 per cent or the total counts. 
As in the spatial analysis the regions marked by green circles in 
Fig.Qare excluded from the spectral analysis. 

For each region, the ancillary response file (ARF) and the re- 
distribution matrix file (RMF), weighted by the X-ray brightness 
in the [0.3-2 keV] energy range, are computed by using the CIAO 
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Figure 2. The surface brightness profile (in units of photon counts) of the 
simulated cluster Cr^i . Open squares represent the values extracted from 
the X-ray analysis, the horizontal bars correspond to the bin sizes. The 
dashed curve is the corresponding best-fit /3-model, with f) = 0.52 and 
r c = 21 kpc (= 0.013 R v i T ). The bottom panel shows the ratio between the 
profile from the X-ray analysis and the best-fit /3-model. 



tools mkwarf and mkrmf . Source spectra are extracted from the 
event file, re-binned and analyzed in the [0.5-6 keV] band. Back- 
ground spectra are extracted from the background event file for the 
same source regions. A thermal model (mekal) absorbed by the 
Galactic column density is fitted to the data by using the \ 2 statis- 
tic in the XSPEC package lArnaudl 19961) . The only free parameters 
are the gas temperature and the normalization, being Galactic ab- 
sorption TVh, redshift z and metallicity Z fixed to the input values 
adopted in the X-MAS run: jV H = 5 x 10 20 cnT 2 , z = 0.175 and 
Z = 0.3Z Q . 

In Table |2| we quote the best-fit values of the spectroscopic 
temperature, T2500. with the corresponding 68% confidence level 
error 02500 ■ These value are also shown as vertical lines in the right 
panels of Fig. Q 



4.3 Deprojection results 

To compute the mass through the equation of the hydrostatic equi- 
librium as described in the next section, we need to recover the 
three-dimensional profiles of the gas temperature and density by 
deprojecting the quantities measured in the X-ray spectral analy- 
sis. W e adopt the deprojection technique presented in Ett ori et alJ 
(2002), that is described briefly here. For each annulus with lumi- 
nosity L r i n g, a single thermal model is fitted to the projected spec- 
trum as described in the previous section, giving the temperature 
T r j ng . By using the geometrical corrections discussed in lKriss et alJ 
( 1983), the photon-weighted, projected-on-the-sky measured quan- 
tities are converted to the values expected in spherical shells by the 
relations: 

1 1/2 



n c = [(Vol T )- 1 #( J B//0.82) 
e=(Vol T )- 1 #L llag , 
T=(Vol r r 1 #(£Wr ring )/e, 



(2) 



where we define the electron density n c in terms of proton density 
n p as n c = n p /0.82, the Emission Integral is EI — J n e n p dV = 
0.82 J n 2 e dV = K x 47rd 2 ng (l + zf x 10 14 , K is the normaliza- 
tion of the mekal model, and the symbol # represents the "ma- 




r/Rvir 

Figure 4. Deprojected density profile of Cr b i\. The red dashed line rep- 
resents the density profile as obtained from the hydrodynamic simulation 
(psim)- The dots are the values extracted from the X-ray analysis (p), the 
vertical bars (having sizes comparable to the dots) are la errors, while the 
horizontal ones correspond to the bin sizes. The bottom panel shows the rel- 
ative differences between the true values of the simulation and the derived 
ones, C = (p sim - p)/p sim . 



trix product" operator. The matrix Vol contains the values of the 
fractions of the spherical volume of a shell seen at each ring. [The 
notation (Vol T ) _1 indicates that the matrix is firstly transposed 
and then inverted]. The outputs are then the measurements of the 
electron density and plasma temperature in each volume shell with 
given inner and outer radii. 

We compare in Fig|3] the deprojected temperature measure, 
T, with the three-dimensional mass-weighted estimate, T mw , re- 
sulting from the simulations. The mass-weighted temperature is 
estimated directly from the hydrodynamic simulations as T mw = 
/ mTdV/ J mdV, where m is the mass of each gas particle. T mw 
is the proper temperature value we sh ould use i n the hydrostatic 
equilibrium equation to derive the mass I M athiesen & Evrardl200lL 
see also Section 4.4). Since the clusters in our sample have an az- 
imuthally quite symmetric thermal structure, we find that their T mw 
profile is not very dissimilar from the T spcc one, even if in the outer 
regions the values obtained in the X-ray analysis are systematically 
lower. A more quantitative comparison can be done by looking at 
the bottom graphs of each panel, where we show a parameter simi- 
lar to the one defined in Section l4~2l A = (T — T mw )/a- Bpcc , where 
now the temperatures are three-dimensional quantities, and the ra- 
tio between the deprojected temperature and the mass-weighted 
temperature: B = (T/T mw ). We find that, while \A\ < 3 up to 
-R500 for all the objects, B indicates that the spectroscopic tem- 
perature is withi n 20 per cent of T mw , wi t h implications on the 
mass estimates jMafhiesen & Evrard^OOlt IVikhlinin et alJl2005L 
see Section 5.1). 

Just as an example, in Fig.|4]we compared to the true profile 
Psim, of the simulated cluster Cueii (dashed line), with the gas 
density obtained from the deprojection technique used. We find a 
good agreement for r/R v - a > 0.1, result that holds also for the 
other clusters. 



Systematics in the X-ray Cluster Mass Estimators 7 




R500 



R500 



R500 



Figure 3. Deprojected temperature profile for Cn ot , Cp er t, C^ergi Cfleii. and Cji e i2- The dashed line represents the profile for the three-dimensional 
mass-weighted temperature T mw as obtained from the hydrodynamic simulation. The open squares are the values extracted from the X-ray analysis, the 
vertical bars are lc errors (cr Bpcc ), while the horizontal ones correspond to the bin sizes. The bottom graphs in each panel show quantities related to the 
differences between the two temperatures: A = (T — T mw )/(T S pec and B = (T/T mw ). The dotted lines indicate A = (—3, 3) and B = (0.8, 1.2). 
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4.4 X-ray estimates of the gravitational mass 

The "true" mass profile, Af s j m (< r), of the simulated objects, ob- 
tained by summing all the masses of the particles inside a sphere 
of radius r, can be now compared to several different X-ray mass 
estimators, M cs t, that we discuss below. 

• Most from the direct application of the hydrostatic equilibrium 
equation: Mhe 

[ hi s method (as discussed in lEttori et alj|2002h makes use of 
the hydrostatic equilibrium (HE) equation under the assumption of 
spherical symmetry to estimate the total gravitating mass: 



A/he(< x) 



kT{x)xR vil fdlnp d\nT 



Gfim 



Rvii f 

p V 



d In x dlnx 



(3) 



where G is the gravitational constant, /i=0.59 is the mean molec- 
ular weight in a.m.u., m p is the proton mass, k is the Boltzmann 
constant, x = r/i? v j r and T and p are the deprojected 3D profiles. 
The equation is applied directly to the deprojected data by estimat- 
ing the derivatives with respect to lnx of ln(p x T) without any 
further smoothing of the measured data points. The HE equation is 
based on simple assumptions: sphericity, static gravitational poten- 
tial and isotropic velocity field. Recent studies have shown how the 
latter hypothe sis is often not s atisfied inside both simulated clus- 
ters (see, e.g.. iRasi a et all2004UK~av et alj20 04|) and o b served ones 
(see, e .g., IZhang et alJ2005tlDupke & Bregmarl2005l) . lRasia et alJ 
(2004) suggested that in order to correctly compute the total mass, 
the gas velocity has to be taken into account, because it is a rele- 
vant component of the total energy equipartition between gas and 
dark matter. The more general model describing the gas dynamical 
equilibrium can be then written as: 



Mhe,v{< x ) 



kT(x)xR vil 
Gpm p 

dlnp 



dlnp dlnT 



G 



dlnx dlnx 
d In a 2 



+ 



dlnx dlnx 



+ 2f3 v (x) 



,(4) 



where /3 V (= 1 — of /2a 2 ) is the gas velocity anisotropy parameter 
and ay and <j t are the gas radial and tangential velocity dispersions, 
respectively. Since some years, velocity fields have been studied 
in sim ulated cluster as bulk motion (see, e.g. Norman & Brvan 
1999) or turbulent motions (see, e.g. lDolag et all2005tlVazza et alJ 
2006). Nevertheless, robust observational measurements of j3 v and, 
more in general, of the gas bulk motions are difficult (see, e.g., 
DuDke & Bregm arl2005h and due to the failure of Suzaku satellite 
spectrometer, we have to wait the new generation of high spec- 
troscopic resolution X-satellites to determine the cluster velocity 
structure. In the present analysis, we evaluate the second part of 
Eq. |4]directly from the simulations. 

• M cs t using (3—model and polytropic temperature profile: 

Assuming that the gas density profile is described by a j3— model 
and that it is related to the temperature through the polytropic rela- 
tion ( T oc p J ~ , with 1^7^ 5/3 ), Eq.|3|can be written as (see, 
e.g., [Henriksen & Mushotz kvl 19861) : 



M 3]7 (< as) = 



kT(x) x fdlnp + dlnT 



Gpm 



-( 

p \ 



d In x dlnx 



3 T o r c x c 
Gpm p (1 + xl) c 

0.757 x 10 14 



/' 



(1+x 2 )" 



Mr. 



(5) 



Table 3. Deviations between the masses calculated by using the different 
mass estimators M cs t (see their description in the text) and the true mass 
M a ; m as directly obtained from the simulated clusters, in units of the es- 
timators errors: |M os t — M s i m |/cr os t. In the case of discrepancies greater 
than 1 Cost we report in parenthesis the value of the percentage difference: 
(|M es t — M s i m |/M cs t) X 100. The values are computed at R2500 and 
extrapolated to R500 for Mg, Mp tl , Mrtm. and Mnfw- 



Cffot 



R2500 

Mhe 11111 

M H E,v 11111 

M Bil > 10 (47) > 10 (49) > 10 (23) > 10 (37) > 10 (36) 

M > 10 (43) > 10 (58) > 10 (24) > 10 (35) > 10 (31) 

M NFW 1 1 1 2(20) 3 (24) 

M RTM 2 (13) 1 1 2(20) 5 (32) 



M 



13,1 



> 10 (52) 



Mp 10 (44) 
Mnfw 1 
Mrtm 3 (56) 



4(45) 
10(51) 

1 

1 



R500 

7 (28) 
10(17) 

1 

1 



9 (42) 7(37) 

410(31) 410(25) 

2 (35) 3(46) 

3 (27) 4(81) 



where x c = r/r c , r c is the core radius (in units of h^Q Mpc), 
q = 1.5/3(7 — 1) + 1 and To is the central temperature in keV. 

• M ea t from an isothermal (3— model: 

Rather than relating the gas density and temperature accordingly 
to the polytropic law, an isothermal gas is here assumed. The total 
mass is then obtained from Eq. |5|by imposing 7 = 1 and To = 
T500 (reported in Table|2j. 

• M e st through analytic mass models: Mnfw and Mrtm 
This method, widely adopted in the reconstruction of the X-ray 

mass profile (see, e.g., lAllen et al . 2002), makes use of functional 
forms of the radial mass distribution obtained from numerical sim- 
ulations. The derived gravitational potential is combined with the 
deprojected gas density profile to recover a temperature profile 
through the numerical inversion of the HE equation. A x 2 distribu- 
tion is then obtained by a grid of values for the two free parameters 
of the mass model (the scale radius r s and the concentration pa- 
rameter c) by estimating the deviation between the numerical tem- 
perature profile and the deprojected one. The errors on the best-fit 
parameters are inferred from the distribution of the \ 2 values. 

As functional f orms of the mass pr ofile, we have used here 
those proposed by iNavarro et alJ 119971) (hereafter NFW) and by 
IRasia etaljj2004l) (hereafter RTM), which read 

M NF w(< x) oc log(l - xcnfw) — XCNF y / j 



Mr 



i(< x) 



{x crtm + 2) 



1/2 



x crtm + 1 



(6) 



respectively, with x = r/7? v i r . 

The values of the corresponding best-fit parameters for our set of 
simulated clusters are reported in Table|5| 



5 RESULTS ON THE MASS PROFILES 

The ratios between the mass profiles reconstructed by using all 
these methods, M es t, and the true mass profiles as extracted by 
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Figure 5. Ratios between the mass profiles derived from the X-ray analysis, M es t , and the true mass profile of the simulated cluster, M a j m for our galaxy 
clusters: Cijot, Cp er t, Gjn erg , Cn e n and Cn e i2- The vertical lines indicate -f?2500- The red asterisks and the green diamonds represent the mass derived 
by assuming the hydrostatic equilibrium, Mjje (Eq.[3j and Mhe v (Eq.|4}, respectively. The solid blue and the thick cyan lines refer to the /3— model, Mg, 
and /3— model plus polytropic relation, Mg _ (Eal5l. respectively. The green and magenta dashed lines show the masses derived by assuming the analytic 
profiles, Mfjpw and Mrxm. respectively (Eq[5|. The measured errors (1 a ) on the mass estimators are represented by vertical error bars or thin solid lines 
in the case of A/nfw an d AferM ■ 
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simulated clusters, M s i m , are shown in Fig. [5] Note that the uncer- 
tainties relative to the X-ray mass measurements as shown in Fig. [5] 
decrease indeed in relation to the dynamic state of the cluster, with 
the largest errors being associated to the first three objects (Cnot, 
Cpert, Cmcvq) and the smallest ones to the most relaxed systems, 
Crcii and Cr^ii- The vertical dashed lines refer to the radii cor- 
responding to i?2500- In Table|3|we quote the differences between 
the true mass and the mass estimates in units of the error at this 
radius and at R500, where Mp, Mp t ^, Mrtm, and Mnfw need to 
be extrapolated. In the case of discrepancies greater than 1 a error 
we report also the percentage difference in parenthesis. 

We summarize here our results and, where needed, discuss the 
main reasons for the observed discrepancies between the estimated 
X-ray mass and the total mass distribution. 

• The mass profile obtained from the HE equation, Mhe, de- 
pendent upon the measured temperature profile (see Fig.[3J- Large 
errors on, and irregular radial distribution of, the temperature val- 
ues induce large scatter on the reconstructed Mhe measurements. 
As clear from the figure, this mass estimator produces large er- 
rors mainly related to the background level. For this reason most of 
these measurements are within la from the expected values, even 
in objects that are still undergone to some (minor) mergers and for 
which the hydrostatic equilibrium might have not been reached yet 
ipHot and Cpert)- If we focus only on the best fit value, we notice 
that, in the case of Cnot, Mhe overestimates the true mass due 
to the steepness of the temperature profile at -R2500 owing to the 
wrong measurement of temperature profile in the last bin. On the 
other hand, it underestimates the true mass in relaxed objects (see 
Fig.EJ. 

The effect of the gas bulk motions are always smaller 10 per 
cent, excluding for the cluster Cp er t- In this case the kinetic energy 
is significant mostly because of the entering merging blob in the 
South and can still contribute between 20 and 30 per cent to the 
total mass measurements. The neglected kinetic pressure term can 
partially compensate for the present best fit deviations in relaxed 
systems. 

• The analytic masses, Mnfw and Mrtm, are reconstructed 
starting from the hydrostatic equilibrium equation, therefore they 
follow the behaviour of the Mhe profile, but they have smaller 
errors. Nevertheless, we notice that the statistical errors are still too 
high to show significant mass discrepancies (Fig.[5}- 

• The /3-model gives the worse mass reconstruction with mass 
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; Rasia et al. 
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Mp n is 5-10 per cent lower than Mp in 
the outer regions (r > R500) owing to the use of a poly tropic law to 
describe the relation between the measured temperature and density 
profiles. Here we note that, in our simulated clusters, the correla- 
tion between In T(r) and In p(r) cannot be reproduced by a single 
linear relation and requires the polytropic index 7 = 1 + lnT/lnp 
to vary between 1.1 and 1.2 moving outwards. The use of an in- 
appropriate assumption on the functional form of the temperature 
profile propagates to the mass measurements making them less ac- 
curate than the estimates done with the isothermal /3-model. The 
main sources of the difference for these models are (i) the poorness 
of the /3-model in describing the density profile, (ii) the incorrect 
assumption of isothermality or above all of the polytropic relation, 
and (iii) the uncertain determination of the /3-model parameters 
(which we will discuss in the next subsection). 

Looking at the behaviour of each single cluster, we notice that, 
even for the dynamically disturbed objects, the central values of 



the total mass can be recovered within 20-30 per cent at -R2500, the 
situation being worse at -R500 • 



5.1 On the NFW concentration and /3 parameter 

The computation of the mass directly depends on the values of the 
best-fit parameters, like the NFW concentration, cnfw, or the /3 
value in the /3— model. The former is evaluated from the mass func- 
tional form that better reproduces the observed temperature profile 
at r > 50 kpc. The inner 50 kpc are also excluded in the spa- 
tial analysis that provides an estimate of /3, to avoid any influence 
from the central cooling region where a too large a mount of cool 
gas co ncentrates in our simulated clusters (see, e.g., iBorgani et alJ 

We investigate here how the measurement of these parame- 
ters depends upon the considered radial range. When we also in- 
clude in the analysis the central region within 50 kpc, the con- 
centration parameter becomes larger, but the reduced x 2 also in- 
creases. The cluster Cmztq is an example: cnfw rises from 8.6 to 
28 and the reduced \ 2 from 0.9 to 2.6. On the contrary, the value 
for /3 does not vary significantly (the most relevant deviation is for 
Cpert, where /3 changes from 0.56 to 0.47). Also in the spatial fit 
of the surface brightness, when including the inner 50 kpc region 
the reduced \ 2 of the fit increases. While this systematic effect can 
be easily handled from an observational point of view, the influ- 
ence of the outer end of the radial range is limited by the field 
of view of the detector and by the background. From a direct fit 
of the data from the simulations, we find that both parameters de- 
pend strongly upon the adopted outer radius, and that, while cnfw 
decreases, /3 increases when larger radii are considered (see also 
[Navarro et al Jl 19951 iBartelmann & SteinmetJfl996t IBorgani et alJ 
2004 . The same trend occurs when we fit an NFW and a /3— model 
to the density profile. For example, we measure cnfw ~ 5.3 by fit- 
ting the simulated density profile of Cm erg out to the virial radius 
(see Table 0, while cnfw ~ 7.0 if we limit the fit to the region 
mapped in the Chandra field-of-view and cnfw ~ 8.6 if we use 
the deprojected temperature profile. 

In summary, the concentration parameter and the /3 value are 
affected by both the method used to derive them and the radial in- 
terval over which they are measured. These two limitations have 
the same effect: the observed cnfw and /3 will be always higher 
and lower, respectively, than the values directly derived by fitting 
the density of simulated clusters up to the virial radius. 



5.2 Effect of the background level 

In the previous analysis, we follow an observational approach to 
recover the mass from Chandra observations. Mainly due to the 
relatively high background level of Chandra , this results in stat- 
ical errors for the temperature profile that, in some cases, are too 
large to be able to show significant mass discrepancies between the 
recovered and true mass in Fig. [5] In prevision of future missions 
that will be able better control and reduce the instrumental noise, 
we reanalyzed the same clusters using Chandra observations for 
which the background accumulation time has been reduced by a 
factor of 100 implying that the effective background has been re- 
duced by a factor of 100. Using these new data set we are able to 
extract cluster signal to larger radii than before (we remind that we 
restricted our analysis to all the annuli with at least 5000 net counts 
and with a flux from the source greater that 25% of the total flux) 
that easily extend beyond i?soo. As before, we find that the all the 
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Table 4. The same of Table |3*1 but derived using a background level 100 
lower that the Chandra value (see text in Sect. 5. 2). 
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best fit values of the estimated X-ray mass underestimate the true 
mass by a remarkable amount: by 20-40 per cent using Mp and by 
< 20 per cent using Mhe or the analytic functional forms. Due 
to reduced statistical error, we can now see that these discrepancies 
are in many cases statistically significant to more than 3cr as clearly 
shown in Table|4| It is worth saying that, in contrast with Table|3|in 
this case the values at -R500 are not extrapolation but measurements. 

It is also worth noting the behaviour of the analytic masses 
in comparison with the case of normal background reproduced 
(Table [5}: being the errors largely reduced the true mass is al- 
ways out of 3ct and it comes to be strongly underestimated. The 
masses recovered by the hydrostatic equilibrium present the same 
characteristic, but for the clusters showing a perturbed dynamical 
state because they are still conserving large error bars. In particu- 
lar, we investigated the reason of the understimation of Mhe and 
we find that half of the total discrepancy is provided by neglect- 
ing the kinetic energy still present as bulk motions of the intra- 
cluster medium. The remaining deviation arises from a systemat- 
ically lower measurement of the representative gas temperature. 
In other words, in these conditions of very low background, the 
deprojected spectral value is systematically lower than the mass- 
weighted estimate from numerical simulations by about 10 per 
cent, suggesting that a proper use of well-defined gas temperature 
values in the hydrostatic equilibrium equation should allow a more 
rigorous measurements of the total mass profile. 



6 SUMMARY AND DISCUSSION 

We have investigated the bias on the X-ray mass estimates with 
mock long Chandra exposures of five galaxy clusters obtained from 
high-resolution hydrodynamic simulations using an observational- 
like approach. The five objects have a spectroscopic-like tempera- 
ture between 2.7 and 11.4 keV, and are at redshift z — 0.175 to fit 
their Rsoo within Chandra ACIS-S3 chip. They are characterized 
by different dynamical states: two objects are perturbed, one had 
a recent major merging, and two are fairly relaxed. These simu- 
lated clusters are processed with our X-ray Map Simulator, X-MAS 
, to produce realistic X-ray event files and images that are analyzed 



with the goal of measuring the total mass profile. To evaluate the 
systematic effects present both in the assumptions of hydrostatic 
equilibrium and of the /3-model and in the observational techniques 
adopted in recovering the mass measurement, we compare the esti- 
mated mass profiles to the ones directly measured in the hydrody- 
namic simulations used as input for our analysis. 

The main results, shown in Fig.|5|and quoted in Tables|3|and 
[4] can be summarized as follows. 

• Due to the relatively high statistical errors mainly connected 
with the high background level of the Chandra observations, the 
HE equation seems to recover, within the errors, the true cluster 
mass profile of the simulated clusters. This partially holds also for 
the analytic mass models NFW and RTM. If however we reduce 
the background level by a factor of 100, we immediately see, to 
a confidence level > 3a, that the reconstructed masses using this 
technique are underestimated by 20% (see Table [4)- In this case, 
we find that the neglected kinetic pressure term can compensate for 
about half of the observed deviations in relaxed systems, while the 
other half is due to the underestimated measurement of the temper- 
ature in comparison to the mass-weighted one. 

• The mass measurements reconstructed via the /3-model are 
the worst among the models considered in the present work since 
they show a systematic underestimate, with typical deviations of 
about 40 per cent at -R2500 and -R500. Mp tJ is 5-10 per cent lower 
than Mp in the outer regions (r > -R500) owing to the use of 
a too simplistic polytropic equation state that relates temperature 
and density profiles. The use of the polytropic functional form of 
the temperature profile makes the mass measurements less accurate 
than that provided from the isothermal /3-model. The other sources 
of biases for the /3-model mass come from (i) its inaccurate de- 
scription of the density profile, (ii) the poor determination of the 
parameters describing the spatially-extended X-ray emission. 

• A systematic effect, that influences the determination of the 
mass through the /3-model mass and the analytic formulae, is also 
produced by the radial interval considered to constrain the value 
of j3 of the /3— model and the concentration parameter in the NFW 
formula. We stress that the observed /3 and cnfw will be always 
lower and higher, respectively, than the values derived directly by 
fitting the density of simulated clu sters out to the virial radiu s (see, 
e.g. th e observational evidence in IVikhlinin et alll99jj : lNeumannl 
l2005t) , 

• We conclude that the mass estimates based on the hydrostatic 
equilibrium equation in combination with the temperature profile 
and those of the analytic fits (NFW or RTM) provide a more robust 
mass estimate than the ones based on the f3— model. 

It is important to remarke that the data analysis performed in 
this paper has been made under ideal conditions. In fact we as- 
sumed to know precisely both the background and the instrument 
response. Furthermore we used very long exposures for all the sim- 
ulated observations reaching a large number of net counts. It is clear 
that the uncertainties in the background and instrument response 
which are inevitably present in real observations (as well as shorter 
exposures) may easily make these discrepancies larger. Moreover 
we remind that the physics adopted in the simulations biases our 
results on the temperature and density profiles more significantly 
when the thermal structure of the X-ray emitting plasma is more 
complex, in the sense that it becomes problematic and ambiguous 
to represent the average properties of a gas which presents a wide 
range in emission measures and temperatures. 

Finally, we notice that these systematic effects on the X-ray 
total mass propagate to the constraints on the cosmological param- 
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eters obtained by using clusters as probes of the matter distribution 
in the Universe. At the present, statistical uncertainties associated 
to the paucity of available high-redshift cluster samples are compa- 
rable to systematic uncertainties arising from biases in cluster mass 
estimates (see, e.g., iRosati et alj|2002|) . as those discussed in this 
paper. As larger samples of distant clusters will be available in the 
near future, it is clear that such systematics will start to dominate 
over statistical errors. In this respect, extensive analyses, like the 
one presented here, will be of crucial relevance to quantify possible 
biases in any procedure of cluster mass estimate. However, the size 
of systematic biases and uncertainties, as calibrated from hydrody- 
namic simulations, depends on the physical processes included in 
the simulations themselves. This highlights the fundamental role 
played by an accurate numerical treatment of the complex physics 
of the ICM to calibrate clusters from hydrodynamic simulations as 
precision tools for cosmology. 
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